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CELLULAR COMMUNICATIONS SYSTEM RECEIVERS 

This invention relates to cellular communications 
system receivers, and is particularly concerned with such 
receivers for use in mobile radio systems, such as TDMA-(time 
5 division multiple access) systems. 

Background 

It is well known that it is necessary in a receiver 
of a cellular mobile radio system to recover each communicated 
signal under varying and challenging conditions. These 

10 conditions include, for example, the presence of multi-path 
signals and fading resulting in low signal-to-noise ratio 
(SNR) , the presence of strong co-channel interfering (CCI) 
signals, and Doppler effects due to relative movement of the 
signal transmitter and receiver, as well as typical constraints 

15 due to factors such as limited channel bandwidth and equipment 
tolerances. In the case of a TDMA system, it is necessary to 
recover the timing and synchronize to the time division 
multiplex (TDM) frames and time slots of a received digital 
communications signal. It is also desirable to provide the 

2 0 receiver with the least possible cost and computational 

complexity . 

International Publication Number WO 97/08867 dated 
March 6, 1997, in the name of Northern Telecom Limited and. 
entitled "Timing Recovery And Frame Synchronization In A 
25 Cellular Communications System", discloses a method of timing 

recovery in which indirect variables of a linear complex vector 
are estimated using a maximum likelihood criterion in order' to 
recover sampling delay and hence the timing and frame 
synchronization of the received signal. 

3 0 There remains a need to provide improved receivers 

for cellular communications systems. 
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Summary of the 'Invention 

According to one aspect of this invention there is 
provided a method of processing samples of a received signal to 
produce a demodulated signal, comprising the steps of: 
representing an impulse response matrix for the received signal 
using a plurality of indirect variables of a linear complex 
vector; synchronizing to the received signal samples in 
dependence upon the indirect variables; tracking the indirect 
variables for successive received signal samples; equalizing 
the successive received signal samples in dependence upon the 
tracked indirect variables; and producing the demodulated 
signal in response to the equalized received signal samples. 

Thus in accordance with this method the indirect 
variables track variations in distortions, such as delay, 
fading, and phase distortions, and their use is extended to 
signal processing steps for producing the demodulated signal, 
thereby facilitating an improved performance of the entire 
receiver of a communications system. 

The received signal may in particular be a signal of 
a TDMA communications system, and the step of synchronizing to 
the received signal samples can comprise matched filtering the 
received signal samples to produce the plurality of indirect 
variables, and determining a maximum of a function of the 
indirect variables to determine synchronization. In a 
particular embodiment of the invention described below, there 
are four indirect variables and said function is a function of 
only two of the indirect variables. 

The step of tracking the indirect variables for 
successive received signal samples can comprise recursively 
filtering initial values of the indirect variables, established 
during the synchronizing step, in dependence upon the 
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successive received signal samples, and can also comprise a 
step of estimating frequency offset in dependence upon the 
successive received signal samples. This enables the tracking 
to be effective over time slots in a TDMA system operating at 
high frequencies, for example 2.4 GHz, despite rapid changes 
due to Doppler effects arising from relative movement between a 
transmitter and a receiver of the received signal . 

The step of equalizing the successive received signal 
samples can comprise adaptively changing an equalizer parameter 
in dependence upon tracking errors for successive received 
signal samples to reduce co-channel interference with the 
received signal . 

The method can advantageously be applied in a dual 
antenna receiver arrangement in which said indirect variables 
are produced and tracked individually in respect of samples of 
a received signal from each of two spaced antennas, received 
signal samples from the two antennas being combined and 
equalized in dependence upon a combination of the indirect 
variables in respect of the two antennas. 

Another aspect of the invention provides apparatus 
for producing a demodulated signal from samples of a received 
signal, comprising: a synchronization unit responsive to the 
received signal samples for producing a linear complex vector 
comprising a plurality of indirect variables having initial 
values corresponding to a synchronized state; a tracking unit 
responsive to the initial values of the indirect variables and 
to the received signal samples to produce tracked values of the 
indirect variables for successive received signal samples; an 
equalizer responsive to the tracked values of the indirect 
variables to equalize successive received signal samples; a 
feedback path from the equalizer to the tracking unit to 
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facilitate producing the tracked values of the indirect 
variables by the tracking unit; and a demodulator responsive to 
the equalized received signal samples to produce a demodulated 
signal . 

The synchronization unit can comprise a plurality of 
finite impulse response filters for matched filtering of the 
received signal samples to produce the plurality of indirect 
variables. The tracking unit can comprise a recursive filter 
for recursively filtering the indirect variables in dependence 
upon the successive received signal samples. 

The apparatus can also include a frequency offset 
estimator coupled to the tracking unit for modifying the 
tracking of the indirect variables in accordance with estimated 
frequency offset in dependence upon the successive received 
signal samples . 

For reducing co-channel interference, the apparatus 
advantageously includes a unit, responsive to tracking errors 
determined by the tracking unit for successive received signal 
samples, for estimating an interference correlation matrix to 
adaptively change a parameter of the equalizer. 

The apparatus can include respective synchronization 
and tracking units for samples of a received signal from each 
of two spaced antennas, the equalizer being responsive to the 
tracked indirect variables for both antennas to combine and 
equalize the received signal samples from the two antennas. 

Brief Description of the Drawings 

The invention will be further understood from the 
following description by way of example with reference to the 
accompanying drawings, in which: 
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Fig. 1 schematically illustrates a dual antenna TDMA 
cellular radio system receiver using indirect variables in 
accordance with an embodiment of the invention; 

Fig. 2 schematically illustrates indirect variable 
synchronization and frame synchronization units of the receiver 
of Fig. 1; 

Fig. 3 schematically illustrates a linear indirect 
variable equalizer of the receiver of Fig. 1; 

Fig. 4 schematically illustrates indirect variable 
tracking, frequency offset estimation, and K<D units of the 
receiver of Fig. 1; 

Fig. 5 schematically illustrates a hard limiter of 
the receiver of Fig. 1; 

Fig. 6 schematically illustrates an optional 
interference correlation matrix estimator of the receiver of 
Fig. 1; and 

Fig. 7 schematically illustrates a single antenna 
TDMA cellular radio system receiver using indirect variables in 
accordance with another embodiment of the invention. 

Detailed Description 

Embodiments of the invention are described below in 
the context of a receiver for use in a TDMA cellular system 
compatible with EIA/TIA document IS-54-B: Cellular System Dual 
Mode Mobile Station- -Base Station Compatibility Standard (Rev. 
B) and later documents referred to as IS-136 and IS-136+, 
minimum performance standards for which are specified in a 
document referred to as IS-138. For brevity, this is referred 
to below simply as the TDMA system and the IS-54 standard. 
However, the principles of the invention are also applicable to 
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other TDMA systems and to other types of communications system 
receiver. 

Such a TDMA system requires a receiver having a high 
performance for operation in various radio channel conditions 
which include a frame duration of 20 ms, a channel bandwidth of 
30 kHz, a multi-path delay spread of up to 42 [is (one TDMA 
symbol) with equal powers of the mult i -path signals, Doppler 
frequency up to 200 Hz, and the presence of up to 3 strong co- 
channel interference signals. The receiver is desired to 
provide the best possible sensitivity and multi-path fading 
reception, despite these conditions, with the least possible 
computation complexity and cost. 

It is well known to enhance reception of radio 
channels subject to fading by using two (or more) spaced 
antennas whose respective receive path signals are combined in 
a desired manner, and the first embodiment of the invention 
described below relates to a dual antenna receiver arrangement. 
However, the invention is also applicable to a single antenna 
receiver, as described later below. 

Despite the use of two antennas, a design of receiver 
for operation in the various conditions outlined above presents 
a significant challenge. As indicated above, one such design 
described in publication WO 97/08867 makes use of indirect 
variables to recover timing and frame synchronization of the 
received signal . 

As described in that publication, each component of 
an impulse response matrix G(x) is approximated by a linear 
combination, plus a constant term, of a pair of functions cpi (x) 
and cp 2 (t). Several examples of function pairs are given. This 
leads to introduction of a variable <D n which is a 3 -dimensional 
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complex vector constituted by the transpose of three indirect 
variables <|> lf n, <h,n, and <t>3,n, which are used in the processes of 
timing recovery and frame synchronization. The recovered 
timing and synchronization parameters are then used in 
conventional manner for deriving the content of the received 
signals . 

The present invention also uses indirect variables, 
but does not merely use them for timing recovery and frame 
synchronization. Rather, the invention recognizes that 
parameters such as the sample timing and frame synchronization 
are only means to the end of recovering the content of the 
received signal, and that these parameters do not necessarily 
provide any value for themselves. Instead, the invention makes 
use of indirect variables substantially throughout the entire 
receiver, and then recovers the content of the received signal 
from the indirect variables at the demodulator. This 
facilitates achieving an improved receiver performance. In 
particular, as described further below, the indirect variables 
can be used in tracking channel changes, Doppler and other 
frequency offsets, providing equalization, and also in reducing 
the adverse effects of co-channel interference. 

Referring to Fig. 1, which illustrates a dual antenna 
TDMA system receiver using indirect variables in accordance 
with an embodiment of the invention, two spaced antennas 10A 
and 10B are coupled to respective receiver front end units 12A 
and 12B, each of which includes a radio frequency receiver, 
down converter, sampler, and analog- to-digital converter of 
known form to provide at its output digital complex signal 
samples y£ and y* , the subscript n denoting the sample number 
and the superscript A or B denoting the antenna. These signal 
samples are supplied to a respective one of two indirect 
variable (IN. VAR.) synchronization units 14A and 14B, to a 
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respective one of two independent variable tracking units 16A 
and 16B, and to a linear independent variable equalizer 18 of 
the receiver of Fig. 1. 

The receiver of Fig. 1 also includes, commonly for 
5 the signals of the two antennas, a frame synchronization unit 
20, a frequency offset indirect variable estimator 22, a hard 
limiter 24, a K<P calculation unit 26, and a demodulator 28. 
Signal connections among these various units of the receiver 
are shown in Fig. 1 and are further described below. In 
10 addition, the receiver can optionally include an interference 
q correlation matrix (ICM) estimation unit 3 0 which is shown with 

~- its connections in dashed lines in Fig. 1. 

jf * 

^ In order to understand the further description below, 

54 it is expedient to consider a mathematical background which 

15 leads to such understanding. This consideration is for a two- 

II 

F| path signal model in an IS -54 TDMA system, in which as is well 

^ known signals are communicated using 7i/4 -shifted DQPSK 

(differential quadrature phase shift keyed) signal symbols in 
non-overlapping time slots each of which comprises data 

2 0 symbols, known synchronization symbols (sync word) , and known 
CDVCC symbols. The known symbols, in particular the sync word, 
facilitate determination and tracking of parameters of the 
signal received during the time slot, these parameters 
including for example carrier phase which can vary during the 

25 time slot. 

With sampling as is usual at twice the clock 
frequency, a discrete observation model for the received signal 
samples has the form: 

Yi = Z u i^4^- - T i - kT - + E ^gfv - *2 - *t - + tu (i) 

k = 0 ^ Z ZJ k = 0 ^ Z AJ 

3 0 where yi is the complex observation sample, i denotes the sample 
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number from 1 to 2N+1 in the observed data sequence, s k are the 
known complex symbols in the sync word of M symbols, 

Ui = ^/rf exp(jv|/o ) are unknown complex amplitude-phase 
multipliers for the different independently fading paths m=l 
and m=2, the samples i on each path m having power P™ and the 
average power of each path being half the average signal power, 
T is the symbol spacing or clock frequency period, t x and x 2 are 
unknown delays of the two paths, g(t) is the impulse response 
of concatenated transmitter and receiver filters, given by: 

, v _ f sinfrt / T) Y cos(a7tt / t) ^ 
9(t) "I 7lt / T - (2at/T) 2 JJ 

where a is the filter roll-off coefficient, and J] ± is a noise 
sequence of complex Gaussian random variables with zero mean, 
variance 2a A , and correlation function 2 a n g (((i-j)T)/2) between 
two random variables T]± and r|j . 

The channel model of Equation (1) can be written in 

the form: 



k = 0 



2 A 2) ~ V 2 



where U™ = UiS M are the amplitude -phase multipliers during the 
sync word and s* = s k (s M )' are transformed symbols of the sync 
word. Assuming that the amplitude -phase multipliers are 
constant during the sync word, then Equation (2) can be 
rewritten in matrix form as : 

Y n = Gfo^uJi + G(X 2 )SU^ + H n (3) 
where Y n = jy^ y 2n ... Y 2n+ 2N-2 y 2 n + 2N-iF ^ a (2N+1) -dimensioned 
observation vector, H n = [^n-x r| 2n ••• *l2n+2N-2 r l2n+2N-iF is a (2N+1) - 
dimensioned vector of correlated noise samples, 
S = (s m )'[s 0 s 1 ... s M _! s M F is an (M+l) -dimensioned vector of known 
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symbols, and G(x) is the impulse response matrix given by: 

g(- t) g(- x - T) ... g(- x - MT> 



G(x) = 



g(- x + T / 2) g(- x - T / 2) 
g(- x + T) g(- x) 

g(- x + NT) g(- x + (n - l)r) 



g(- x + T / 2 - MT) 
g(- x - (M - l)T) 

g(- x + (N - M)r) 



a 

m 



m 



According to the IS-54 standard, 

T ± € (- L pr T / 2; L pr T / 2) i = 1, 2 |x x - X 2 | < T 

5 where L pr is the number of sample spacings in an uncertainty 
range of delay, abbreviated below to L. With this uncertainty 
range determined to be L=2 (path delays within two sampling 
intervals) , then to ensure that the sync word symbols are all 
within an observation set the above impulse response matrix 
10 must be increased by two initial rows and two final rows, so 
that with N=M it becomes a matrix with 2M+5 rows and M+l 
columns with components Gij=g ( -x+ (i-l-L) T/2 - ( j -1) T) where i is a 
row index from 1 to 2M+5 and j is a column index from 1 to M+l. 

In publication WO 97/08867 an approximation 
15 g(x + iT / 2) = a 1/i (p 1 (x) + a 2 f i<p 2 ( T ) + a 3,i i = ... - 2, - 1, 0, 1, 2, ... or, in 
matrix form, g(x) = A 1 cp 1 (x) + A 2 cp 2 (x) + A 3 i- s used for a single path 

channel to approximate the impulse response matrix G (x) , where 
ai f i, a 2 ,i, and a 3/ i are approximation coefficients and Ai, A 2 , and 
A 3 are approximation matrices with these coefficients. This 

2 0 approximation is reasonable for values of the variable x within 
the interval [-T/2; T/2], but is not sufficient for two paths 
for which, with frame synchronization established, the 
uncertainty range is still [-3T/2; 3T/2] . A good accuracy with 
this greater interval has been found with four terms, i.e.: 

25 g(x + iT / 2) = a l i cp 1 (x) + a 2 i cp 2 (x) 4- a 3 ^cp 3 (x) + a 4 i (p 4 (x) i = ... - 2, - 1, 0, 1, 2, ... 

or G(x) = A 1 cp 1 (x) + A 2 cp 2 (x) + A 3 cp 3 (x) + A 4 cp 4 (x) (4) 
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Various approximation functions can be used to 
provide a desired approximation accuracy, and the invention is 
not limited to any particular set of approximation functions. 
As one example, the approximation functions may be: 
cp-^x) = sin (tit / 2T) cp 2 (x) = cos (tit / 2T) 

<p 3 ( T ) = sin (7CT / 4T) cp 4 (t) = cos (nx / 4T) 

As another example, the approximation functions may be: 

<PiC0 = ^2 CO = §ft> 

(p 3 (x) = -1 . 2 2g(T) + g(T / 2) (p 4 (x) = -§(t) + 1 . 41g(T / 2) 

where g(x) is the Hilbert transform of the function g(t) . These 
approximation functions provide a signal -to -approximation noise 
ratio of about 3 0 dB in a range of t from -T to T. As can be 
appreciated, a preferred set of approximation functions, and 
the number of functions in the set, depends on the desired 
approximation accuracy (signal-to-approximation noise ratio) 
arid directly affects the resulting complexity of implementing 
the approximation functions in the receiver. 

Because the difference in delays between the two 
paths is not more than one symbol spacing interval, Equation 
(3) above can be rewritten, using the approximation functions, 
in the form: 

Y n = A^cp^K + cPi(x 2 >Jn) + A^cp^K + cp 2 (T 2 )j2) ^ 
+A 3 S(( P3 (t 1 )uJ 1 + cp 3 (T 2 >4) + A^cp^K + cp 4 (T 2 )u£) + H n 

If a 4 -dimensioned vector <D n of indirect variables is 
defined by: O n s (^n J 2( n' ^n' Knf 

where <|> ifn = + q>i(* 2 >£ i = 1,2,3,4 (6) 

are the indirect variables, then the model of Equation (5) can 
be expressed in the form: 

Y n = B<D n + H n (7) 
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where B = [(a x s) (a 2 s) (a 3 s) (a 4 s)] is a known matrix because Ai to A 4 
comprise fixed coefficients and S is the known sync word. 

In this context, in the receiver of Fig. 1 the 
indirect variable synchronization units 14A and 14B serve to 
produce initial values 0>£ and <D* respectively of the indirect 
variable vector O n for synchronization, and the indirect 
variable tracking units 16A and 16B serve, in conjunction with 
the frequency offset indirect variable estimation unit 22, to 
track the indirect variable vector <D n throughout a time slot to 
produce tracked values o£ and <D* respectively. The linear 
indirect variable equalizer 18 comprises a Kalman filter which 
is controlled by the tracked values <I>£ and <3>* o f the indirect 
variable vector to combine and recursively filter the received 
signal samples y£ and y* , thereby producing a received and 
equalized signal vector S 0 , n . This vector is limited by the 
hard limiter 24, from the output of which the K<D calculating 
unit produces a feedback control signal Ko, n for the tracking 
units and the demodulator 2 8 produces a demodulated signal on 
an output line 32 of the receiver. The various units of the 
receiver are further described below. 

Fig. 2 illustrates one form of the indirect variable 
synchronization unit 14A (the unit 14B is similar) and one form 
of the frame synchronization unit 20, for producing the initial 
value Oq for synchronization. The unit 14A comprises four 
finite impulse response (FIR) filters (FIR-1 to FIR-4) 40 which 
are supplied with the received signal samples y£ , a calculation 
unit 42, a delay line 44 comprising delay elements each 
providing a delay of one symbol spacing interval T, and a 
selector 46. The frame synchronization unit 20 comprises a 
combiner 48, a delay line 50, and a maximum detector 52. 
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20 



In order to simplify matrix inversion, the model of 
Equation (7) is divided into even and odd sample sets so that 
the model can be expressed as: 

v odd _ „ ^ , rrodd y even „ <T) + H even 

Y n B odd <P n + H n Y n ^even^n ^ n n 

5 and, because B is a known matrix, the indirect variable vector 
can be initially determined by a matched filtering represented 

^odd Lt „ V 1 ^ v odd ^even _ Lt r y even 
by: 0£ aa = ^ dd B odd j B odd Y n O n - V B even B even/ B even% 

in which the first three terms on the right-hand side of each 
equation can be pre-calculated and stored. It is observed that 
10 here and below the various signal processing operations produce 
results which are estimates rather than the precise values of 
the respective signals. In Fig. 2, the FIR filters 40 perform 
this matched filtering function, the filters FIR-1 and FIR-3 
corresponding to the first and third rows respectively of the 

15 matrix (b^B^^B^Y^ and the filters FIR-2 and FIR-4 

corresponding to the second and fourth rows respectively of the 

B even B even ) B even Y n 



,even 
3 even *n 



Consequently, the outputs of the FIR filters 40, 
which are supplied to inputs of the delay line 44, constitute 
the indirect variable vector <D n in accordance with the above 
model, but its synchronization, i.e. the value n which provides 
a reference timing point, is not yet determined. This is 
determined as described below using, for simplicity, only the 
first two of the approximation functions described above, the 
25 outputs of only the filters FIR-1 and FIR-2 being supplied to 
the calculation unit 42, which calculates and produces at its 
output a value modO n . The calculation carried out by the unit 
42 is dependent upon the particular approximation functions 
which are used as described above. For the functions <p x (x) = g(x) 
and cp 2 (x) = g(x) referred to above, for example, 



30 
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[<P 2 frn)i>2,n + cPi^n^nFk^n^.n + <Plfrn))>i,n] 



2 

where x n = — real 
n 



( fx. \ 



a tan 



4>2,n 
♦l,n. 



and the calculation unit 42 



calculates mod<D n accordingly. 

The values modO^ and thus produced for the two 

paths A and B are combined by the signal combiner 4 8 in the 
frame synchronization unit 2 0 to provide for frame alignment of 
the signals from the two antennas, and the output of the signal 
combiner is supplied to the delay line 50 having a length 
corresponding to the observation window of the received signal 
samples. The maximum detector 52 determines the value n 
corresponding to a maximum one of the outputs of the delay line 
50, thereby determining synchronization, and supplies this to 
the selector 46 in the unit 14A. The selector 46 selects from 
the delay line 44 the corresponding indirect variable vector <D n 
and supplies this to its output as the initial indirect 
variable vector <Dq . 

One form of the linear indirect variable equalizer 18 
is illustrated in Fig. 3 and comprises a Kalman filter 
including signal combiners 60 and 62, multipliers 64 and 66, 
and a delay element 68, and an arrangement for determining 
filter parameters including calculation units 70, 72, 74, and 
76 and a delay element 78. The arrangement and operation of 
the equalizer will be further understood from the following 
description. 

With prediction of estimates to an n-th step at the 
step n-L, the Kalman filter in Fig. 3 can be seen to provide 
its output S 0/ n in accordance with: 
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S 0in = A n S n _! + K S(n (y n - G{O n / ^Xa^^)) 
where K s , n is a Kalman filter gain given by: 



(8) 



K S/ n = fen I n-X®n / n-jfefan / n-Jfen / n-A®n / n-J] + 

V S ,n/n-l = A n V s#n _X + V 5#n O) 

V S , n = Vs^/n.i - % Jv s#n , / n _ L )J 

where the terms of these equations can be understood from the 
following description. 

Extending a single antenna representation for the 
case of two antennas and hence two sampled signals, 



10 



Yn = 




s n + 



Tin- 
Tin 



where the last matrix represents equivalent noise including 
indirect variable errors. This noise has the covariance 

" 1 g(T / 2) 
g(T / 2) 1 _ 



matrix ; 



°2x2 



u 2x2 



15 where V d 



V ri,4> = B ( S n/n-lK,n/n-L B ( S n/n-l) 
"*,n/n-L = V d>,n-L + Q<D L i s the covariance matrix of the 
indirect variable prediction error for L+l steps, Qo is the 
covariance matrix of exciting noise of the indirect variable 



vector, B(s) = 



s t aJ 

S T A* 



, G(O n ) = 



^nAi 



, A x and A 2 are 4 x (2m+l) - 



20 



dimensioned approximation coefficient matrices as described 
above, and m is a number of adjacent symbols taken into 
account. A vector condition for TDMA symbols can be described 
by the equation S n = A^^ + £ n where A n is a shift matrix (or a 
transition matrix during the CDVCC) and ^ n is noise with 
covariance matrix V^ n = 2Q^ n which is a zero matrix during the 
2 5 CDVCC, thus when n is not in the CDVCC: 
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A n = 



0 
0 

0 
0 



1 

0 

0 
0 



0 

1 

0 
0 



0 
0 

1 

0 



Qn = 



0 
0 
0 



0 
0 



and when n relates to a known CDVCC symbol w n : 

0 10 ... 0 
0 0 1 



0 
0 



0 
0 



1 



Qn = 



0 0 
0 0 
0 0 



Thus referring again to Fig. 3, the calculation unit 
70 produces G(O n /n-L), used in the first and third of Equations 
(9) and in Equation (8) , from the indirect variable vectors 
and tracked as described below, and the calculation units 

76, 74, and 72 produce the values of respectively the first, 
second, and third of the Equations (9) . V^ n supplied to the 
unit 74 is predetermined as indicated above, and V n/I1 supplied 
to the unit 76 can also be fixed and predetermined or, as later 
described below, can be adapt ively changed. The Kalman filter 
gain K s , n is used by the multiplier 66 to produce, with the 
other elements 60, 62, 64, and 68 of the Kalman filter, the 
equalizer output vector S 0/n from the signal samples y£ and y* 
in accordance with Equation (8) . 

Referring now to Fig. 4, one form of the indirect 
variable tracking unit 16A is illustrated, the unit 16B being 
similar and connections to it also being indicated in Fig. 4, 
together with one form of the frequency offset indirect 
variable estimation unit 22 and of the K<D calculation unit 26. 
To some extent these units make use of similar calculations so 
that they are closely inter-related and they are accordingly 
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described together below. As illustrated in Fig. 4, the 
indirect variable tracking unit 16A comprises signal combiners 
80 and 82, multipliers 84, 86, 88, and 90, and a delay element 
92. Although not shown separately in Fig. 4, the unit 16A is 

5 also supplied with the initial value <E>^ of the indirect 

variable vector O n for synchronization by the unit 14A as 
described above. The frequency offset indirect variable 
estimation unit 22 comprises signal combiners 94, 96, and 98, 
multipliers 100 and 102, a delay element 104, and calculation 
q 10 units 106, 108, and 110. The K<D calculation unit 26 comprises 

In 



o 
in 

a 



signal combiners 112 and 114, a delay element 116, and 
calculation units 118, 120, 122, and 124. 



jjj The effect of frequency offset in conjunction with 



indirect variable error £, can be expressed by: 
15 0£ = dF n O^ + # n 

< = dF^jU + 4f, n (10) 

dF n = + 5dF,n 

where dF n » exp(j27i^ f T) and f Q f is frequency offset. With an 
observation model of the form: 
2 0 y£ = B(s n >D£ + t£ y* = B(s n )0^ + t,* 

then if dF n _i is estimated very accurately, substituting the 
estimates for the actual variable in Equation (10) produces a 
model which does not depend upon this variable: 

25 <D* = <DjU + 4*,n (10) 

and results in a filtering algorithm: 

< = ^n-l<-l + K*,nkn ~ B(S n >I>£_ldFn-l) (11) 
where = V 0 n / n.^sjv^ 
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V" 1 , = (B(s n K,n/n-lB(S n j # + 2^^]"' (12) 

V <P,n = V 0,n/n-l " K O,n[ V 0>,n/ n-l B ( S n) ] 
V <D,n/n-l = V 0,n-1 + Q^,<P + v dF,n-l 

It can be appreciated here that the same Kalman gain matrix Ko, n 
is used for signals from the two antennas, thereby simplifying 
the filtering algorithm and making its complexity independent 
of the number of antennas . 

It can be seen from Fig. 4 that the elements 80 to 92 
of the tracking unit 16A are arranged to implement Equation 
(11) for the signal samples y£ , the parameter dF n -i being 
supplied from the output of the delay element 104 and the 
parameter Ko,n being supplied from the output of the unit 124, 
to produce the tracked indirect variable vector 0£ at the 
output of the signal combiner 82. The four Equations (12) are 
implemented in Fig. 4 respectively by the calculation units 
124, 122, and 120 and the elements 112 to 116. 

For frequency offset estimation, the following 
approximated model is derived from the above: 

dF n = dF n _! + ^ dF/n 

Yn = B(s n >D^ 1 dF n . 1 + v£ 

y* = B^^^dF^ + v* 
where v£ and are equivalent observation noise which take 
into account estimation errors of the indirect variable vectors 
On_x and and have the same covariance matrix defined by: 

V V/n = B(s n )v 0fn/n _ 1 B(s n )' + 2a^R n 

Then an indirect variable filtering algorithm for dF n 
can be written in the form: 
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dF n = dF n _ x + K^ F#n ^ - B(s n >D^ 1 dF n . 1 )+ K* FfTL (yl - B(s n K_ 1 dF n _ 1 ) (13) 
where K£ F , n = v^.tPn-i B ( s n)' V'n 

K* F ,n = v dF/n O^_ 1 B(s n yv- 1 n (14) 



V dF/ n = 1 / 



l/v d F f n/n-l + + -<&n-l B(S n J V"* B^)©*^ 

v dF,n / n-1 = v dF,n-l + Q^,dF 



Although Equations (14) appear to be very complex, 
they can be implemented with low complexity because a large 
number of steps are already otherwise performed. For example, 

the inverse matrix is used for calculation of Ko,n, the 

differences in Equation (13) are the same as those in Equation 
(11) , and multiplications such as B(s n )<t>n_i have been used for 
difference calculations . 

It can be seen from Fig, 4 that the four Equations 
(14) are implemented by the calculation units 108, 110, and 106 
and the signal combiner 94 respectively. The multipliers 100 
and 102, signal combiners 96 and 98, and delay element 104 
implement Equation (13) . 

Referring to Fig. 5, the hard limiter 24 is supplied 
with the vector S©, n from the output of the equalizer 18, and 
supplies the elements s n+m to s n - m of this vector via respective 
ones of 2m+l stages 130 to produce hard limited elements which, 
with delayed versions thereof produced by delay elements 132, 
constitute the elements of the vector S n that is supplied to the 
KO calculation unit 26, and in particular to the unit 118 as 
shown in Fig. 4. In addition, the hard limiter 24 provides an 
output for the hard limited version of the element s n _ m to the 
demodulator 28, which operates in a well-known manner for 
demodulating the 7i/4-shifted DQPSK signal. Each of the stages 
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13 0, as represented in Fig. 5 for two such stages, provides 
signal phase rotation, hard limiting, and derotation following 
the 7i/4 -shifted DQPSK modulation rules. 

In order to optimize parameters for operation of the 
receiver as described above, particular values can be selected. 
For example, the number 2m+l of simultaneously estimated 
symbols in the equalizer 18 can be selected as being 5 with 
m=2, and with L=2 as already indicated the number of symbols 
m+l+L used in the tracking units 16A and 16B is also 5. The 
integration interval for the synchronization units can be 8T, 
and in Fig. 4 Q^o can be (diag ( 10" 4 * [2 1 0.5 0.5]), and Q^, dF can 
be 2.5*10~ 5 . In addition, a fixed value of SNR, for example 17 
dB, can be used for synthesizing the above algorithms, as the 
actual signal-to-noise and interference ratio may be unknown. 

The receiver as described above is intended to 
provide an optimum performance in the presence of noise. 
However, in the presence of co- channel interference, the 
performance of the receiver can be degraded. In order to 
reduce or avoid such degradation, the receiver can also include 
the interference correlation matrix (I CM) estimation unit 30 
shown in dashed lines in Fig. 1. This provides an adaptive 
control of the matrix V^n which is supplied to the linear 
indirect variable equalizer 18, so that interference 
cancellation is also achieved by the operation of the equalizer 
as described above. 

An analysis can be carried out in a similar manner to 
that described above in relation to the operation of the 
tracking units, but in respect of the information symbols Si in 
the TDMA time slot rather than the synchronization symbols, 
from which it can be determined that the noise covariance 
matrix V n#n provides interference cancellation based on 
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differences of estimated values sf = yf - G(of ^ and 
8? = yf - 0(0?^ which are already determined (Equation (8) 
above, using a slightly different notation) in the operation of 
the Kalman filter as described above. As illustrated in Fig. 
1, these differences are supplied from the tracking units 16A 
and 16B to the 1CM estimation unit 30 to enable estimation of 
the matrix V n , n . This matrix is determined by: 



V n = 2R^ = 2 



2 

a A R rj r AB R ti 
r AB R r| a B R r| 



where a| and are unknown variances (real variables) , r AB is 
an unknown correlation coefficient (complex variable) , and is 
the known covariance matrix already specified above. 

Fig. 6 illustrates one form of the ICM estimation 
unit 30, in which the unknowns , , and r^ are each 
averaged over a desired number NA of samples by respective 
delay lines 140 and summing units 142, the outputs of which are 
supplied to a calculation unit 144 to determine the covariance 
matrix V n in accordance with the above equation. The unit 3 0 
also includes multipliers 146, 148, 150, 152, and 154, and 
transpose units 156 and 158, which serve to produce the 
elements of the matrix as inputs to the delay lines 140 for 
averaging. Thus it can be seen that this adaptive operation of 
the receiver adds very little complexity to the receiver, but 
can substantially improve the performance of the receiver in 
the presence of co-channel interference. 

Although a dual antenna embodiment of the invention 
has been described above in detail, it should be appreciated 
that the invention is not limited in this respect, and it may 
also be applied to a single antenna receiver as illustrated in 
Fig. 7. Thus as shown in Fig. 7, a single antenna 160 is 
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coupled via a receiver front end unit 162 whose output signal 
samples y n are supplied to an indirect variable synchronization 
unit 164, an indirect variable tracking unit 166, and a linear 
indirect variable equalizer 168. The synchronization unit 164 
5 provides an initial synchronization vector O 0 to the tracking 
unit 166, no frame alignment being required because in this 
case there is only one signal path. The tracking unit 166 
provides a tracked vector <D n , and this is equalized by the 
equalizer 168 to produce a resulting signal for demodulation by 
10 a demodulator 170. It can be appreciated that the units 164, 
□ 166, and 168 in the receiver of Fig. 1 can use similar 

3 principles to those described above in order to provide 

J ] improved single-antenna receiver performance through the use of 

jy. 

^ indirect variables for all of the signal processing in the 

tl 15 receiver prior to the demodulator. 

^ It can be appreciated that although as described 

ri 

=? J 

above the frequency offset estimation unit 22 is provided as is 
preferred to compensate for frequency offsets, which may be due 
to local oscillator frequency variations and, especially, due 
20 to Doppler effects, in other embodiments of the invention this 
unit can be omitted. 

In addition, although the description above refers 
to, and the drawings illustrate, particular units such as 
calculation units, signal combiners, multipliers, delay 
25 elements, and so on, it should be appreciated that in practice 
the functions of all of these units can conveniently be carried 
out by one or more digital signal processors or application- 
specific integrated circuits. 

Thus although particular embodiments of the invention 
3 0 have been described above, it can be appreciated that these and 
numerous other modifications, variations, and adaptations may 
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be made without departing from the scope of the invention as 
defined in the claims. 
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